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We will re-derive E.S. formulation for triangular shape charge distribution. 
The triangular shape charge distribution is represented by: 
3a* 1 
(=- irl) for |r| < 1/a 
0 for |r| > 1/a 
This equation represents a spherically symmetric triangular shape charge distribution 


whose magnitude is zero when the distance from the atomic center is more than 1/a. The 
Ewald coefficient parameter is the only parameter that effects the width of the 
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distribution. The magnitude for |r| < 1/a is given as see (2= |r| ). The reason why the 
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triangular distribution is represented in this form comes from the definition of the direct 
sum component. At any point where |r| > 1/a, the net direct sum potential due to the 
point charge and the charge distribution is zero. Therefore, when computing the force or 
potential at a point, we can neglect, without error, all charges at a distance greater than 
r=1/a. 


To verify the result; potential component is computed by: 
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For the Gaussian the integral becomes: 


Since the triangular distribution is defined in two different regions (|r| < 1/a and |r| > 
1/a), we need to derive the potential for these two regions separately. For |r| < 1/a 
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To get the force for the direct sum in this region we simply take the partial derivation of 
the direct sum potential: 
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For |r| > 1/a: 
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And the direct sum force in this region is: 
1 1 
qe (2, ¢-)) = 
ATE, r 7 
And the reciprocal sum potential can be found by: 


1/a 4 . 

dB 3a* 1 Sin[2zkr] 

= 4 2 (= = ) ae 
Or(r) ATE I ( a Na)” 2nkr a 


dp 3a?Sin [| (—knCos [| + aSin [=] 


ATE k*n4 


Comparing the Gaussian and Triangular Charge Distributions 


The preference between the Gaussian and the triangular should depend on the 
acceptable error and computation time. 


The minimum error is possibly achieved by the Gaussian distribution. In the existence of 
a parallel reciprocal sum, if the total consumed time can be decreased to the time 
consumed by the direct sum then we can expect the triangular to perform better than the 
Gaussian at small and large cutoffs (in time) and to perform better in the means of both 
error and time at small cutoffs. 


2/2 


